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Abstract 
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■ We present a chiral-invariant meson-baryon Lagrangian which describes the 

interactions of the baryon octet with the lowest-mass meson nonets. The 
nonlinear realization of the chiral symmetry generates pair-meson interac- 
tion vertices. The corresponding pair- meson coupling constants can all be 
expressed in terms of the meson-nucleon-nucleon pseudovector, scalar, and 
vector coupling constants, and their corresponding F/(F + D) ratios, and for 
which empirical estimates are given. We show that it is possible to construct 
an NN potential of reasonable quality satisfying these theoretical and empir- 
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^ , ical constraints. 
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I. INTRODUCTION 

The construction of AW potentials has a long history, and a large number of different 
models have appeared in the literature. After a comparison of these AW potentials by 
confronting them first with the pp scattering data [0J and later also with the np scattering 
data ||, it was found that all existing models do a rather poor job (some worse than others) 
in properly describing these data. This is rather disturbing, since most of these models 
are used in many-body calculations and pp bremsstrahlung calculations to draw certain 
conclusions on three-body forces, relativistic effects, off-shell effects, etc. And so one can 
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ask the question whether it is at all valid to conclude anything from a many-body calculation 
if the applied NN potential cannot even adequately describe the two-nucleon data. 

This uncertainty has been recently cleared up with the construction of a number of high- 
quality NN potential models [0,0. With high-quality we mean that these models describe the 
NN scattering data with the almost optimal xV-^data ~ 1- Although based on very different 
functional forms, these potentials give remarkably similar results in many-body f)|-[7| and pp 
bremsstrahlung H calculations. One 'disadvantage' of these new potential models, however, 
is that they are largely phenomenological and only explicitly include the one-pion exchange. 
Our next task, therefore, is to try and construct a high-quality potential model which is less 
phenomenological in that it explicitly contains the heavier-meson exchanges (with coupling 
constants which can be compared with empirical data) and which, possibly, is consistent 
with the symmetries of QCD. 

Recently, Ordonez, Ray, and van Kolck || presented a nucleon-nucleon (NN) potential 
based on an effective chiral Lagrangian of pions, nucleons, and A isobars. Using 26 free 
parameters, the agreement with the experimental scattering data was found to be satisfac- 
tory up to lab energies of about 100 MeV. An extension to higher energies and a further 
improvement in the description of the data would require an expansion to higher orders 
in chiral perturbation theory, making the model much more complicated and introducing 
many new parameters. Hence, the authors conclude that it is not practical for potential 
models derived from effective chiral Lagrangians to compete with more phenomenological 
approaches. 

In this paper we would like to advocate an alternative approach. We want to investigate 
whether it is possible to construct a potential model which not only gives a satisfactory de- 
scription of the scattering data up to ~350 MeV using a limited number of free parameters, 
but which also retains the salient features of chiral symmetry. For that purpose we here 
do not integrate out all mesons other than the pion [ 1Q| , but rather adopt the successful 
approach used in one-boson-exchange models and keep all lowest-lying mesons with masses 
lower than 1 GeV, say. The AW potential model is then obtained by evaluating the standard 
one-boson-exchange contributions involving these mesons, but now including the contribu- 
tions of the box and crossed-box two-meson diagrams [[II]] and of the pair-meson diagrams 



where at least one of the nucleon lines contains a pair- meson vertex | 12]j . The potential 
contributions are calculated up to order 1/M 2 in the nucleon mass as is customary in the 
conventional one-boson-exchange approaches, and so we want to stress that here we are not 
doing chiral perturbation theory. Here we only use chiral symmetry to generate the interac- 
tion Lagrangians and to find constraints for the associated coupling constants. However, as 
we will see, the pair-meson diagrams arise as a direct consequence of chiral symmetry. The 
box and crossed-box diagrams then also have to be included because they are of the same 
order in the number of exchanged mesons as the pair-meson diagrams. The pair-meson 
interactions can be viewed |K| as the result of integrating out the heavy-meson (masses 
larger than 1 GeV, say) and resonance (e.g., A, N*, Y*) degrees of freedom in the two- 
meson-exchange processes. Also, according to "duality" fllBf , the resonance contributions 
to the various meson-nucleon amplitudes can be described approximately by heavy-meson 
exchanges. Treating the heavy-meson propagators as constants, which should be adequate 
at low energies, then leads directly to pair-meson exchanges. 

In Refs. |TT|,P~2"[ we already showed that the inclusion of the two-meson (box, crossed-box, 
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and pair) contributions provides a substantial improvement in the description of the scat- 
tering data as compared to a potential containing only the standard one-boson exchanges. 
Although originally [0,0 the meson-pair interaction Lagrangians were taken to be purely 
phenomenological, we later found that the pair-meson coupling constants could all be fixed 



using experimental input and chiral-symmetry constraints. In particular, in Ref. [12] the 
estimates for the pair-meson coupling constants were based on the linear a model |16| . In 
order to appreciate this result, it should be realized that, by fixing the pair- meson coupling 
constants in this way, this improvement could be obtained without the introduction of any 
new parameters. 

This remarkable result encourages us to go beyond the AW model and to investigate 
whether a similar approach will also be fruitful in the construction of an extended hyperon- 
nucleon (YN) potential. An important motivation for the development of an extended YN 
model is provided by the study of hypernuclei using one-boson-exchange models Hl7| , |l8| . The 
problem with the construction of a YN potential is that there are only a few experimental 
data available, and these data are rather old and not very accurate. At present it is very 
difficult (if not impossible) to determine all the free parameters from the scattering data 
alone. A reduction in the number of free parameters is obtained by first fixing the parameters 
which also play a role in NN scattering (and which are easier to determine, since there are 
many accurate NN scattering data), and then using SU(3) symmetry to fix the coupling 
constants in the YN potential. This approach has been used successfully in the various 
hard-core [|I| and soft-core p(J YN one-boson-exchange models of the Nijmegen group. 



The Julich models even impose an SU(6) symmetry. As a further reduction in the 

number of parameters, it would be convenient if we were able also to fix all (or at least most) 
single-meson coupling constants at their empirical values. 

As a first step, in this paper we will construct a chiral-invariant Lagrangian for the 
meson-baryon sector. Rather than extending the linear a model (the model we used in |L2| ) 
to the strange sector, we here find it more convenient to explore the nonlinear realization of 
the spontaneously broken chiral symmetry. This allows us to express all pair-meson coupling 
constants in terms of the single-meson vertex parameters. Using symmetry arguments and 
further experimental data we give estimates for these single-meson vertex parameters. We 
then have two options for constructing a potential model: (i) we can either strictly impose 
all the constraints on the coupling constants and investigate whether it is at all possible 
to impose chiral symmetry in a potential model or, (ii) we can relax at least some of the 
constraints and, hence, optimize the description of the scattering data. As an application we 
will here only explore the first option and show that a fully constrained NN model indeed 
allows for a very satisfactory description of the data. This is a remarkable result because, 
until now, in potential models it has never been possible to fix more than only a few coupling 
constants at their empirical values. The second option reflects that the imposed symmetries 
need not be exactly true, a fact which can be useful when fine tuning the model to the 
scattering data. This option will be left for the future. 

The outline of the paper is as follows. In Sec. we list the baryons and mesons of the 
model. Section |T| gives a review of the linear a model ]16[ and the nonlinear realization 



of chiral symmetry as introduced by Weinberg ||23|| . We briefly indicate how the model can 



be made to agree with experimental data. In Sec. [IV] we extend the nonlinear realization 
to SU(3) to obtain a chiral-invariant meson-baryon Lagrangian. This SU(3) Lagrangian 
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describes the meson-baryon interactions of the scalar, pseudoscalar, vector, and axial- vector 
meson nonets with the baryon octet. Section ^ then lists estimates for the various single- 
meson vertex parameters (singlet and octet coupling constants, mixing angles, and Fj (F+D) 
ratios), using theoretical and experimental input. The coupling constants for the double- 
meson vertices can all be expressed in terms of these single-meson vertex parameters; this 
will be discussed in Sec. [Vl|. Unfortunately, we have not been able to construct a potential 
model that satisfies all these constraints. Therefore, in Sec. |V1 we briefly indicate how 
the interaction Lagrangian can be extended by introducing new free parameters without 
abandoning the empirical constraints on the single-meson coupling constants. This allows 
for much more flexibility in actually imposing these constraints in a baryon-baryon potential 
model. As an application, in Sec. [VI I| we show that with this extension it is indeed possible 
to construct an AW potential that satisfies the constraints of chiral symmetry and which 
gives a very satisfactory description of the AW scattering data up to 350 MeV. 



II. MESON AND BARYON CONTENT 

In the following, the baryons are the members of the SU(3) octet: A(940), A(1115), 
£(1192), and S(1318). The mesons consist of the standard pseudoscalar nonet (ir, rj, rj', K) 
and vector nonet (p, u, <p, K*). The physical isoscalar mesons in these nonets come about 
as an admixture of the pure octet and pure singlet isoscalar states. That is, the physical r\ 
and rj' are admixtures of the pure octet rjg with the pure singlet r/o, and the physical lo and 
<f> are admixtures of the pure octet cug with the pure singlet ojq. 

In the SU(3) xSU(3) chiral- invariant model, the extension of the global chiral symmetry 
to a local chiral symmetry introduces two octets of gauge fields, which can be re-expressed in 
terms of one octet of vector and one octet of axial- vector fields (see below). It seems natural 
to identify the octet of vector gauge fields with the octet of vector mesons. The singlet (which 
is not a gauge field) is then added to complete the nonet. Similarly, the octet of axial- vector 
fields can be identified with the octet of axial- vector mesons. Adding the singlet and mixing 
the pure octet and pure singlet isoscalar states is then assumed to give the physical nonet of 
axial- vector mesons [ax(1260), /i(1285), /i(1420), and ^(1270)]. In this scenario, the axial- 
vector mesons are a necessary ingredient in establishing the local symmetry, but one can 
argue that their role in low-energy baryon-baryon potential models is likely to be negligible. 
Their masses are well above 1 GeV, and so the corresponding short-range potentials are 
probably already effectively included by form factors at the meson-baryon vertices. For 
alternative approaches to include vector and axial-vector mesons in a chiral-invariant way, 
see Refs. @J2J. 

The existence of the ++ scalar nonet (octet plus singlet) is still controversial. The 
Particle Data Group |26| lists an isovector state ao(980) and an isoscalar state / (980), 



while there appears to be evidence for a strange isodoublet of scalar mesons, denoted by 
k(887) f27|. These would make up the octet. There is also evidence for a broad isoscalar 
++ (760) state ||28|| , which could be the scalar singlet. Alternatively, in analogy with the 
pseudoscalar and vector nonets, it is likely that also the /o(980) and e(760) are admixtures 
of the pure octet and pure singlet isoscalar states. The baryon-baryon potential due to 
the exchange of a broad meson can be approximated by the sum of two potentials where 
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each potential is due to the exchange of a stable (effective) meson | 29| . Due to the large 
width of the £(760), the low-mass pole in this two-pole approximation is rather small (~500 
MeV), and hence would represent a possible candidate for the low-mass a meson which is a 
necessary ingredient in baryon-baryon potential models, and which is the ingredient of the 
linear a model. It is beyond the scope of this paper to discuss whether or not a nonet of 
scalar mesons with masses below 1 GeV really exists, how its existence can be explained or 
disclaimed in a quark model, or whether a scalar one-boson exchange is nothing else but an 
effective representation of correlated two-pion exchange. Here we only mention that already 
some time ago Jaffe [31J presented a quark-bag model calculation of q 2 q 2 mesons where 
the members and decay modes of the lowest multiplet fit in remarkably well with what is 
observed for the scalar nonet discussed above. More recently, Tornqvist |H) has calculated 
the properties of a distorted qq nonet in a unitarized quark model, which also identifies the 
above scalar nonet reasonably well (except for the k(887) state). Finally, after an absence 
of almost 20 years, the latest edition of the Particle Data Group now tentatively lists a 
very broad resonance as /o (400-1200) with a width of 600-1000 MeV. Here, we will stick to 
e(760) and assume a width of T e « 800 MeV. 



III. SU(2) CHIRAL SYMMETRY 

Although the SU(2) case has been discussed extensively in the literature, we believe it 
is still instructive to first review the SU(2) case in some detail before we turn to the SU(3) 
case. In the SU(2) case the results can be easily expressed in terms of the physical fields, 
which makes it easier to see what is going on. In the SU(3) case the expressions are much 
more involved. 

The classic example for a model with chiral SU(2) symmetry is the linear a model flTB" 



The linear a model contains an isotriplet of pseudoscalars, n, and an isosinglet scalar, a, 
which can be grouped into 

£ = cr + IT -7T, (3.1) 

and which transforms under global SU(2) L xSU(2)/j as 

E -> LERl (3.2) 

The nucleon field if) has left and right components, ipi,R — 1(1 T Ts)^ transforming as 

ip L -> L^ L , ^ R -> R4> R . (3.3) 

Given these transformation properties, we can construct the chiral-invariant Lagrangian 
£ = C^p + £2 + £/, where 

C s = iTr(^St^ S ) _ I M 2 Tr(£t£) _ lA 2 Tr(£t S ) 2 , 

£ I = -g$ L ^ R + ^ R Z^ L ), (3.4) 

and n, A, and g free parameters. Obviously, the ground state cannot have the full 
SU(2) L xSU(2)£. symmetry, since in that case all hadrons (nucleon, pion, sigma in this 
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model) would have partners of equal mass but with opposite parity, which is not the case 
in the real world. Hence, the chiral symmetry has to be spontaneously broken down to the 
vectorial subgroup SU(2)y. This can be achieved by adding a term linear in the a field, 
choosing the ground state as (a) = f , and shifting the scalar field to a = s + f . We then 
find the familiar Lagrangian 



C = ^(i^d^ — M)ip — gijj(s + ijsT-7r)ip 

+§(d"«d„s - m 2 s s 2 ) + K^tt-^tt - m 2 7r 2 



m; 



2/o 



S{S 2 + 7T 2 ) 



m; 



*/o 2 



{s 2 + * 2 )\ 



(3.5) 



where the nucleon mass is given by M = gfo, and the scalar and pseudoscalar masses are 
given by m 2 = /i 2 + 3A 2 /o an d = /2 2 + A 2 /q, respectively. The PCAC (partially-conserved 
axial-vector current) condition for the axial-vector Noether current, d^A^ = mlf^TT, deter- 
mines the relation between (a) = f and f v = 92.4 ± 0.3 MeV |2B|, the pion decay constant. 
As we will see below, the introduction of vector and axial-vector fields enforces a renormal- 
ization of the pion field, and so we cannot make the identification fo = (a true identity 
in the absence of vector and axial- vector gauge fields). 



Weinberg has shown |32[] that the linear a model has the major disadvantage that it hides 
the fact that soft pions are emitted in clusters by derivative couplings from external lines, 
and so it is more convenient to transform the nonderivative ipi'jsT -TZip and an 2 interactions 
away. For that purpose, Weinberg defines a nonlinear transformation of the nucleon fields, 
given by |3| 

(3.6) 



^ = (l+^)-l/2(l-, 75T .£)iV, 

with £ chosen such that 

ip(M + gs + ig^T-Tr)^ = N(M + gs')N. 
It will be convenient to write this transformation in the form 



N R 
N L 



l + iT-j 

(1+ £2)1/2- 



in which case we find 



E' = f + s' = u ] T,u ] . 



(3.7) 



(3.8) 



(3.9) 



Solving Eq. ( |3.7|) gives 



7T 



f0 + S + VOb + ^ + TT 2 



= — —TT 



s' = [(fo + s) 2 + n 



211/2 



fo 



2/o 

(tfE) 1/2 -/o, 



(3.10) 



where the chiral rotation vector £ is proportional to a new pion field, it' '. We should mention 
that the condition ( |3.7| ) also removes the air 2 interaction, and the new Lagrangian only 
contains three-point and four-point interactions for the s' scalar field. 
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In this nonlinear realization of the spontaneously broken chiral symmetry, the new fields 
u(£) of the coset space SU(2) L xSU(2)^/SU(2)y transform as [[33 



u -> LuH ] = HuR\ (3.11a) 
H = VLu 2 WRu^ = VRuWLu, (3.11b) 



where the equality in Eq. ( p. 11a ) is due to parity. The left and right components of the new 
doublet nucleon field N both transform with the same SU(2) matrix H, and so the nucleon 
part of the Lagrangian can now be written as 

C N = W(i'fD IJi — M)N — gNNs - gjJt'f'fupN, (3. 12) 

where we defined = + iT^, and 



4/ 2 i + e 

1 r-d^ir' 



Ufl = -f (u%u - ud^) = — — (3.13) 



The transformation rule for the connection T^, 

r M — > HVfj,H^ — iHdfj,H\ (3.14) 

ensures the invariance of the kinetic-energy term. The transformation rule for the N field, 
iV — > HN, means that NN is already an invariant in itself, and so the mass M can be 
treated as a free parameter. Similarly, the transformation rule for u^, 



■ft 



)1 ' " u.^ 



ff« M fft, (3.15) 



allows for the introduction of the free parameter in Eq. (|3.12 ). Choosing it to be 



gA = 1.2601 ± 0.0025 [26], which is the value for the weak interaction axial-vector cou- 
pling constant, and substituting for ir', the pion is seen to couple to the nucleon via the 
pseudovector coupling with strength gA/2fn- This relation between the weak axial- vector 
and the pion pseudovector coupling constants is known as the Goldberger-Treiman rela- 



tion [[34]]. Although this relation only exactly holds in the chiral limit, experimentally it is 
found to hold within about 2%. Note that in the ip representation of Eq. (|3.5|) , this coupling 
is found to be gA — 1, which is off by 25%. 

The next step is to add the isotriplets of vector (p) and axial- vector (ai) fields. One 
way to do this is to extend the global chiral symmetry to a local one and to define the 
corresponding left- and right-handed gauge fields as 



In — \ T ' V — \ T ' (pfi + a /i) j 

r» = = |t-(p m - a M ), (3.16) 



with their field strength tensors 

Ip, = dpl v - dJu + ig v [l,t, Q , (3. 17) 
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and similarly for r^ v . Their transformation properties are 



9v 

r M -> Rr^B) - —Rd^R\ r> -> ifr^, (3.18) 
9v 

and so the proper field strength tensors for the nonlinear transformation are uH^u and 
ur,,„U- Invariance under local chiral transformations implies the extensions 



u\d^ + ig v lp)u + u(d^ + ig v r^)u ] 

u ] {d^ + igyl^u - u{dp + igvr^u^ . (3.19) 

The locally chiral-invariant Lagrangian can now be written as [35] 

C N = NiiYdf, - M)N - g s NNs' - \g v N^TN-p'^ - ^XN^rN-a!^, (3.20) 

where g s , gy, and A are free parameters, and where we defined the field abbreviations 

2 2 1 

a i = ^ ~ TTIP2 * x ^ x ® + —TT~^ ^ + x Pn)- (3-21) 

These field combinations are seen to give rise to pair (itit, np, irax) vertices, and other 
higher-order multiple-meson vertices. 

As an example of how this model can be made to agree with the actual experimental 
results, we next give some attention to the meson sector. The axial-vector field can be given 
a mass by adding to the Lagrangian a chiral-invariant term proportional to Tr(it M w M ); this 
will also generate the kinetic-energy term for the pion field. The vector field receives its 
mass from a term proportional to Tr(/ M /^ + r M r At ), which is only invariant under global chiral 
transformations. This term also contributes to the mass of the axial-vector fields. However, 
in order to get an acceptable agreement with experiment (not only with respect to the 
empirical masses, but also with respect to some of the empirical meson coupling constants 
and decay widths), we have to extend the Lagrangian even further. A very general form in 
the context of the linear a model is given by Ko and Rudaz [36|. Translating their results 
in terms of the fields of the nonlinear transformation ( |3.8|) , the masses for the vector and 
axial- vector mesons are obtained from the Lagrangian 

C m = / 2 Tr(n%) + \ rr%Ti(l% + r%) - cg 2 v (fo + s') 2 Tr(Fu 2 r^ 2 ). (3.22) 

We find that there is a mixing between the it and a M fields, which can be removed by 
making a field redefinition for the a M field, a M = — h(d^Tz + <?y7rx p), and then choosing 
h appropriately. Also, the kinetic-energy term for the pion field is no longer of the canonical 
form, and so we have to do a wave function renormalization: n = Z~ l l 2 Tt r . The final result 
reads (for details, see Ref. |36| ) 
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2 2 2 f2 

m P = m o ~ c 9vfoi 
m l, = m l + ( C + I)0v7i 



ai — ""o ^ ■^yy./O) 
^ = l-^/oV<, (3-23) 

and the PCAC condition gives f n = Z^ 2 f = 92.4 MeV. The parameter c is needed to 
get a simultaneous agreement for the p mass (m p = 770 MeV), the a\ mass (m ai = 1.23 
GeV), and the pNN coupling constant (gy = 5.04 from p° — > e + e~); all data from Ref. ||26|| . 
Substituting f Q in the equation for leads to the two solutions (Z n , c) = (0.173, —0.131) and 
(0.827, 1.25). Similarly, the numerical value for the pirn coupling constant (g p7T7T = 6.05 from 
p — > 7r + 7r _ [F2BJ) can be enforced by adding yet another term to the Lagrangian, introducing 



a new parameter k 6 (see Ref. [0). Having fixed both c and Kq, some of the other meson 
properties such as the root-mean-square pion radius and a\ decay widths also come out very 
favorably p6[| . 



Although a full discussion of the meson sector is outside the scope of this paper, it 
is important to note that it is indeed possible to construct a chiral-invariant Lagrangian 
(sometimes at the expense of introducing new parameters and small chiral-symmetry vi- 
olating terms) which does a reasonably good job in describing a variety of experimental 
results. Furthermore, the introduction of the axial-vector field as a gauge boson enforces a 
renormalization of the pion field, and so we explicitly need part of the Lagrangian for the 
meson sector to define this renormalization constant. 



IV. SU(3) CHIRAL SYMMETRY 



The extension to SU(3) is easily obtained by replacing the Pauli isospin matrices, Tj, by 
the Gell-Mann matrices, A a . The two-component nucleon field is replaced by a 3x3 traceless 
baryon-octet matrix where the left- and right-handed components transform as 



(4.1) 



The meson content of the model still consists of the isosinglet scalar, a, but the pseudoscalar 
triplet is extended to an octet by including the 7]$ and the four K mesons. They are 
collectively denoted by n a . These nine mesons are grouped into 



S = a + iXaTTa, 



1,... 



(4.2) 



which still transforms according to Eq. (|3.2|) , but where L and R are now elements of SU(3). 
Translating the definition ( |3.8j ) to the SU(3) case, we write 



Br 
B L 



u(£ a ) = exp[zA a £ a ] = exp 



2/o 



(4.3) 



where we used a more convenient representation for tt(£ a ), the elements of the coset space 
SU(3)z / xSU(3)/j/SU(3)y. We will identify the primed fields with the octet of physical pseu- 
doscalar fields, which is given by 
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y/2 y/6 



71 



rO 



V 



7T 



1T 



^0 



if 

v/6 y 



(4.4) 



and where the fields transform according to the SU(3) analogue of Eq. ( |3.11a| ). Following 
the phase convention of Ref. [37], the new octet of baryon fields reads 

\ 



B 



/ E° A 



A 

V2 VE 



V 
n 

2A 



(4.5) 



which transforms as 



B -»■ 



(4.6) 



The slightly different transformation property of the baryon octet matrix requires that the 
SU(3) analogue to the nucleon derivative, D^N = (d^ + iT^N in Eq. ( |3.12|) , reads 



D lt B = d lt B + i\T lt ,B). 



(4.7) 



Due to the transformation properties ( [4.1D of the original octet fields, we cannot simply 
copy the interaction Lagrangian of Eq. ( |3.4|) to the SU(3) case. To ensure invariance under 
chiral SU(3)lxSU(3)r, the S field would then have to appear quadratically; i.e., the inter- 
action Lagrangian has to be of the form Tr(ty L J}ty rY^ + ^/^S^lS). Alternatively, we note 
that the field combinations uY$u and vtHifl both transform according to Eq. ( f4.6|) . Defining 
the combination 



we thus find that the most simple interaction Lagrangian is given by 

Cj = -g. A Tr(B~x+B) - g s , 2 Tr(BB X +), 



(4.8) 



(4.9) 



where g Sj i and g s ^ are arbitrary constants. In principle, we also could have included a term 
of the form gpTr^B^X-B), but that would have re-introduced a nonderivative pseudoscalar 
interaction Q which we wanted to get rid of in the first place. 



1 This in contrast to the SU(2) case, where we can choose u(£j) such that \- = 0, and so all 
nonderivative pseudoscalar interactions are transformed away. The difference is due to the fact 
that (Tj 7Tj ) (tj TTj ) is proportional to H2, whereas (A a 7r a )(A;,7rf ) ) is not proportional to I3. 
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Unfortunately, the interaction Lagrangian ( |4.9|) is not good enough if we want it to 
generate the empirical baryon masses, since it gives the same mass, M = (g S) i + g s ,i)fo-, for 
all baryons in the baryon octet. This problem can be solved by adding an octet of scalar 
fields, A a cr a , where the isoscalar octet member is given a nonvanishing vacuum expectation 



(4.10) 



value. We can then extend £ by writing (Ao = y 2/3U3) 
E = F + X s + X a s a + iA a 7r a , (a = 



with 7r a the original pseudoscalar fields and F the vacuum expectation value of the scalar- 
nonet fields, i.e. 



(4.11) 





( fx 








F = 





fl 







V 





h 



with 



fl 



(°"o) + \A(o- 8 ), h 



(a ) - 2J\(a 6 ) 



(4.12) 



In the Appendix we show that now it is still possible to find a transformation «(£), which 
transforms away the octet of original pseudoscalar fields while leaving the vacuum expecta- 
tion matrix F invariant, and where the £ a fields can be identified with an octet of new pseu- 
doscalar fields. As before, this transformation generates the combinations x± °f Eq. ( |4.8| ). 
Both combinations contain the original pseudoscalar fields in a complicated way. However, 
X+ behaves like a set of scalar fields, and so we can simply define these new fields to be the 
physical scalar fields and drop any reference to the original scalar fields. Clearly, to zeroth 
order in the new pseudoscalar fields, the old and new scalar fields are the same. Dropping 
primes, the nonet of new scalar fields is given by 



X+ = F + A s + A a s a , 
where sq now denotes the new scalar singlet and the octet is given by 



(4.13) 



—A, 

V2 a a 



y/2 V6 



\ 



V 



Q _j_ J^8_ 

^0 



K 



2£g 

\/6 J 



(4.14) 



The pure octet isoscalar s$ and pure singlet s are mixed to give the physical / (980) 
and e(760). Similarly, — i\- can be identified as a new isosinglet pseudoscalar field, not 
present before. Because it transforms as H{—ix-)H\ it can be formally added to the octet 
pseudoscalar matrix, which completes the nonet. It must be remembered, however, that 
group transformations such as given for instance in Eq. ( |3.11a|) and which involve different 
transformation matrices, are only valid for the (traceless) octet matrices. 

As discussed for the SU(2) case at the end of the previous section, also in the SU(3) case 
the global chiral symmetry can be easily extended to a local chiral symmetry. The required 



11 



gauge fields are now given by two octets of combinations of vector and axial-vector fields. 
The vector octet is given by 



V2 V& 



\ 



p~ 

K* 



p^ J\ 
\/2 V6 



\ 



K 



*o 



Ve ) 



(4.15) 



and a similar matrix for the axial- vector octet. Adding a mass term of the form (|3.22|) breaks 
the local symmetry and introduces a mixing between the pseudoscalar and axial- vector fields, 
which requires a redefinition of the axial-vector fields. This in turn means that the kinetic- 
energy term for the pseudoscalar fields is no longer of the canonical form, which is taken 
care of by renormalizing the pseudoscalar fields. Unfortunately, then it is not possible to get 
a satisfactory agreement for both the vector and axial-vector masses simultaneously. But 
also this problem can be solved. One way is to introduce yet another term |36| proportional 
to Tr^^Er^E^), which renormalizes the vector and axial- vector fields. Another option is 
to modify the term proportional to Wq by inserting the combinations EE^ and E^E. The 
advantage of the latter extension is that it does not involve a new parameter, while it does 
allow for a very satisfactory description of both vector and axial- vector masses. Therefore, 
a convenient spin-1 mass Lagrangian reads 



(4.16) 



Using this Lagrangian, we find that the set mo = 653 MeV, c = —0.131, f± = 222 MeV, and 
f 2 = 290 MeV gives very satisfactory results for the p, ai, K*, and Ki masses, as well as for 
the pion and kaon decay constants. Given these parameters, the renormalization constants 
are Z n = 0.173, = 0.225, and Z V8 = 0.236. Note that (Z n ,c) is the same as in the 
SU(2) case, and so the good phenomenology in the p-ai sector |36|] is still valid. To make 
a further connection with the real physical world, we next also introduce isosinglet vector 
and axial-vector fields, which complete the experimentally observed nonets. The isosinglet 
fields are taken to be SU(3) invariants. This means that we can formally include them on 
the diagonal in the respective octet matrix representations. 

Finally, to be complete, we should also list the kinetic-energy and mass terms for the 
scalar and pseudoscalar fields. The kinetic-energy term for the scalar fields is given by 
Tr (D' J 'x+D fl x+)i with the covariant derivative defined as in Eq. ( f4.7|) to ensure the chi- 
ral invariance. The kinetic-energy term for the pseudoscalar fields is already contained in 



Tr (u^u^). The simplest possible mass term for the scalar fields is of the form 
A = f?Tr( X+ A s ) - c 2 / 1 2 Tr(xU + ) - c 4 Tr( X U+) 2 , 



(4.17) 



where the first term with the diagonal matrix A s = diag(x, x, y) breaks the chiral symmetry. 
This term has to be included in order to remove the terms linear in the Sq and s§ fields as 
generated by the last two terms. This condition fixes (x, y) in terms of (c 2 ,c 4 ). Fitting to 
the scalar masses as discussed in Sec. |, we find (02, c 4 ) = (7.62, —0.46). 
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To generate the octet pseudoscalar masses, we again need a term that breaks the chi- 
ral symmetry. The singlet pseudoscalar mass is generated by Tr(x_X-)i which is chiral 
invariant. Hence, an appropriate Lagrangian for the pseudoscalar meson masses looks like 



i m 2 

4 "So 



(4.18) 



with the diagonal matrix A p = diag(x', x', y'), where x' = m^/Z n and y' = 2m 2 K / Z^— m%/Z n . 
Note that with this choice the quadratic Gell-Mann-Okubo mass formula, 3m^ g + m% — 
4m 2 K = 0, is approximately satisfied. The mass m Vo is chosen such that, with the proper 
mixing angle, we get the empirical 77 and rj' masses. 

The main achievement of the field transformations as discussed above is that we have 
constructed various 3x3 matrices of the general form $ = (l/y/2)\ c <f> c , where c = 0, ... ,8. 
These matrices contain scalar, pseudoscalar, vector, and axial-vector meson fields, and (ex- 
cept for the vector mesons) they all transform in the same way as the baryon fields. This 
allows us to define the following chiral-invariant combinations 



BB<$> 
BB$> 
~BB<$? 



D 



= Tr(B$B) -Tr(BB<f>), 

= Tr(£$£) + Tr(££$) - § Tr(55)Tr($), 

= Tr(573)Tr($). 



(4.19) 



For a traceless matrix $ (i.e., only octets of mesons), the 5-type coupling [BB&] g van- 
ishes, and the F- and D-type couplings can be written as Tr(_B[$,_B]) and Tr(_B{$, B}), 
respectively; a notation often encountered in the literature. Defining a baryon-baryon-meson 
octet coupling constant g oct and a baryon-baryon-singlet coupling constant g sm , a general 
interaction Lagrangian which satisfies chiral symmetry can now be written in the form 



d = -g oct V2{a [bB<$> 



a) 



BB<S> 



D } 



BB§ 



s' 



(4.20) 



where a is known as the F j [F + D) ratio, and the square-root factors are introduced for 
later convenience. The various field matrices are given by 



y/2 



[F + A c s c ] , 



V2g v 



"7m 



u 1 



V2g v 



u ] (d^ + \gvK{(f + A"- hD^) c ) v 
+u (> + \g v \ c {(f -A» + hD»n) c ) 
V (&* + {gvK{ff + A»- hD^) c ) u 
-u (> + y v X c (ff -A» + hD^n) c ) 7/ 



(4.21a) 



(4.21b) 



(4.21c) 



where D^(Xn) = <9 M (A7r) — ±g v [(\n), (Ap M )], and h chosen such that the mixing between the 
axial- vector and pseudoscalar fields in the meson sector vanishes. Note that the pseudovector 
coupling of the pseudoscalar fields is already included in the axial- vector field matrix $ ax . 
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In addition to the electric coupling 7^/?^ [where p 1 is a shorthand for the fields appearing 
in Eq. (|4.21b[) 1, it is also possible |35j to include a chiral-invariant magnetic coupling (T^ u p', 
where p' is the field strength tensor for the p' field combination. This is due to the fact 
that p' transforms according to Eq. (|3.14j) , and so we can define a chiral-invariant field 
strength tensor p' , as in Eq. ( |3. 17| ) . The transformation ( |3.14| ) also imposes the constraint 
that the chiral- variant D-type coupling [BBQ vc ]d should vanish, i.e., the electric a v = 1. 



This represents the so-called universality condition proposed by Sakurai PBfl . Hence, the 
assumption that the p meson couples universally to the isospin current in this model is a 
direct consequence of chiral SU(3) symmetry. On the other hand, the magnetic a™ is still a 
free parameter. 

Until now, it appears that we have not gained much by imposing the SU(3)xSU(3) 
chiral symmetry: the form of the interaction Lagrangian ( f4.20| ) can be written down imme- 
diately by assuming only an SU(3) symmetry, and has been known for a long time (see, e.g., 
Ref. |3£|). However, the important difference is that the 3x3 matrices $ are not just simple 
representations of the scalar, vector, or axial-vector meson fields, but they contain the pseu- 
doscalar fields in a nonlinear way as well. This means that the chiral Lagrangian contains all 
kinds of multiple-meson (pair, triple, etc.) interactions not envisaged before. The coupling 
constants for these multiple-meson interactions can all be expressed in terms of the single- 
meson interaction coupling constants. This will be the subject of Sec. |V]]. Furthermore, 
the pseudoscalar fields are renormalized due to the introduction of the axial-vector fields. 
Hence, already in leading order the chiral Lagrangian gives interactions between baryons 
and pseudoscalar mesons that are slightly different from what is obtained with the standard 
(i.e., nonchiral) Lagrangian. Also, the coupling constants of the pseudoscalar mesons are 
directly related to those of the axial- vector mesons. Finally, the imposed symmetry allows 
us to fit all the octet meson masses with only two parameters for each octet. 



V. SINGLE-MESON VERTICES 

We will first look at the single- meson baryon-baryon coupling constants, i.e., the interac- 
tion terms linear in the meson fields. Let us drop for a moment the Lorentz character of the 
interaction vertices (IL4 in spinor space for scalar mesons, 757^^ for pseudoscalar mesons, 7^ 
and a^d" for vector mesons, and 757^ for axial-vector mesons) and take as an example the 
nonet of pseudoscalar mesons. The derivative (pseudovector-coupled) pseudoscalar-meson 
interaction Lagrangian to lowest order in the fields is then of the form 

£pv — £p V + £pv^> (5-1) 
where the S"-type coupling in Eq. ( |4.2(J| ) gives the singlet interaction Lagrangian 

m n £$ = -fNN V0 (NN)ri - f SEvo (EE) Vo - fex Vo (S ■ S)r/ - f A A Vo (AA)r/ , (5.2) 
with the (derivative) pseudovector coupling constants 

fNN V0 = /bE7) = /s£»70 = fAAvo = Ip™, (5-3) 
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and where we introduced the charged-pion mass as a scaling mass to make the pseudovector 
coupling / dimensionless. The interaction Lagrangian for the meson octet is obtained by 
evaluating the F- and D-type couplings in Eq. ( |4.20| ), and can be written as 

m^f = -/^(iVriV)-7r-/ H ^(SrS)-7r-/ A& (AS+m)-7r + ^ EEw (SxS)-7r 
-fNN V8 (NN)r] 8 - f B s m ( EE )Vs - f AAvs (AA)r] 8 - / ES??8 (S-S)r/ 8 
-Iank [(NK)A + A(KN)] - f 3AK [(EK C )A + A(K C E 

-fe NK h-(KrN) + (NtK)-'Z] - / 3EK \E-(K c tE) + (HriQ-E . (5.4) 



Here we introduced the doublets 



N 



K 



and S and 7r are isovectors with phases chosen p9| such that 





S-7T = S H 



7T 



]°7T 



7T 



(5.5) 



(5.6) 



The octet coupling constants are given by the following expressions (/ = f° 

-^/(l + 2a), 

-^/(4«-l), 

/(I — 2a), 



fNNTr 



■f, 



2a), 



/as^ = 7= /(I - a) 



EStt 



2/a, 



fNN m 
/AAr/8 : 



73 /(4a -1), 
-^/(l + 2a), 

-^/(1-a), 



Iank 
fsAK 

fllNK 

/hex 



(5.7) 



Similar relations (without the scaling mass m n ) are found for the coupling constants of the 
scalar, vector, and axial-vector mesons. Note, however, that for the pseudoscalar mesons the 
relations ( f5.3| ) and ( |5.7| ) should be slightly modified due to the difference in renormalization 
factors Z n , Z K , Z rv and (see Sec. |TV| ). 

For each type of meson there are only four parameters: the singlet coupling, the octet 
coupling, the F/(F + D) ratio, and the mixing angle to generate the physical isoscalar 
mesons from the pure octet and singlet fields. In most cases we can impose theoretical 
and experimental constraints on these parameters. This will be discussed next. As already 
mentioned before, the axial-vector mesons are very heavy, and hence are not expected to 
play an important role in low-energy potential models, but we will still discuss them here 
for reasons of completeness. 



A. Scalar mesons 

Because the existence of a nonet of scalar mesons with masses below 1 GeV is still contro- 
versial, the constraints on the SU(3) parameters solely depend on the particular theoretical 
model one wants to use to describe the scalar mesons. As a matter of fact, it is not at all 
clear whether it is indeed valid to impose an SU(3) symmetry. However, in order to limit the 
number of free parameters, we will here assume the standard SU(3) relations and assume 
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that the physical e(760) and /o(980) mesons are admixtures of the pure octet as and pure 
singlet ctq fields, in terms of the scalar mixing angle 6$, 



\f ) = sin 6g |<7g) + cos 6» 5 |cr ), 

\e) = cos 6s |cr 8 ) — sin 8 S \a ). (5.8) 

A possible value for the mixing angle is obtained by assuming that the scalar mesons are 
q 2 q 2 states pOfl , and that the e(760) does not contain any strange quarks (hence, its low 
mass). This implies an ideal mixing angle 9s = 35.3°. 

A value for as can be determined by assuming that the baryon masses are generated 
by the nonvanishing vacuum expectation value of the <7o and erg scalar fields. Defining the 
vacuum expectation values fi and / 2 as in Eq. ( |4.12| ), we have 



v^cro) = (2/ a + / 2 ), V3(a 8 ) = (/j - / a ), (5.9) 
or, assuming the ideal mixing 9s = 35.3°, 

fi = (/o(980)), / 2 = -V2<e(760)>. (5.10) 



From the interaction Lagrangians (|5.2| ) and ( |5.4|) for the scalar fields we find the following 
relations for the baryon masses, 

M N = M - Kf(4a s - l)(.f 2 - fi), 

M A = M -fC t («5-l)(/ 2 -/i), 
M E = Mo + |^(a s -l)(/ 2 -/i), 

M H = M + ^ s ° c ct (2a 5 + l)(/ a - A), (5.11) 

with M = C n (2/i + ^/v^ = |(M A + M s ). According to the relations given above, 
these masses satisfy the equality 2M^ + 2M= = 3M A + Ms, which is indeed approximately 
true experimentally (4516 MeV versus 4537 MeV). Solving for as, we find cts = 1.42 and 
C ct (/ 2 - fi) = 136.5 MeV. 



Finally, using the estimates for f\ and / 2 as given in Sec. [TV], we find gf™ = 3.8 and 
<? s °c ct = 2.0. 



B. Pseudoscalar mesons 



The mixing angle 9ps for the pseudoscalar mesons is defined by 

\r)) = cos 9 p S \r] s ) - sm9 PS \r] ), 

\t]') = sin 9 p S \r) 8 ) + cos 9 PS |^o) • (5. 12) 

The linear and quadratic Gell-Mann-Okubo mass formulas give 9ps ~ —23° and 9ps ~ 
— 10°, respectively. The current experimental evidence, however, seems to favor |4(| w 
-20°. 

The axial- vector current F/D ratio, obtained from the Cabibbo theory of semileptonic 
decays of baryons, gives |TJ 0.575±0.0165, or a PS = 0.365 ± 0.007. 
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The ttNN pseudovector coupling constant is obtained from the Goldberger-Treiman rela- 
tion f34] , which gives = /mvtt = fl'A(0)m 7r /2/ 7r = 0.952 ± 0.003. It can also be extracted 



from multienergy partial-wave analyses of pp, np, and pp scattering data [|2[ , which yields 



for the coupling constant at the pion pole f% NlT /^ = 0.0745 ±0.0006. This latter result has 
to be extrapolated to t — before comparing it with the Goldberger-Treiman result. 

Finally, we can also give an estimate for the singlet pseudovector coupling constant. The 
estimate is based on the fact that we can rewrite the 7/ 8 and t]q in a nonstrange-strange basis, 
rather than the standard {u, d, s} quark basis, and then assume that the purely strange- 
quark state does not couple to the nucleon. To be specific, in the standard quark basis 



1 778) = J\\uu + dd — 2ss), 

\rj ) = J h\uu + dd + ss) , (5.13) 



and so the purely strange-quark state can be expressed as \S) = \ss) = (\rjo) — v^l^s)) / 
This state does not couple to the nucleons provided that /atat^ = Inn^ or > equivalently, 



C. Vector mesons 

The mixing angle 9y for the vector mesons is defined by 

\u) = sin 9y \u>s) + cos 9 V \u> ), 

\(f>) = cos 6*y \cus) — sin^y |u;o). (5-14) 

Ideal mixing (i.e., the meson is a pure ss state) gives 9y = 35.3°, which is very close to the 
experimental value (6y ~ 35°) and the values found for the linear (8y ~ 36°) and quadratic 
(6 V ps 39°) Gell-Mann-Okubo mass formulas |26[ . 



As mentioned earlier, the transformation property of the connection requires that 
a v = 1 for the electric coupling of the vector mesons, which is the universality assump- 
tion §§. 

The pNN coupling constant is obtained from the electromagnetic decay p° — > e + e~. 
The vector-meson dominance (VMD) hypothesis assumes that this decay proceeds via the 
photon and that the p°-7 coupling is proportional to the pNN coupling f38fl . The decay 
width m T(p° -> e + e~) = 6.77 ± 0.32 keV then gives c/° c ct = g NNp = 2.52 ± 0.06. Note the 
factor of 2 difference between the definition of in Eq. ( |4 . 2 0| ) and gy in Eq. ( |4.21b| ). 



If we assume that the meson is a pure ss state, which does not couple to the nu- 
cleons (i.e., an ideal mixing angle 9y = 35.3°), then we find for the singlet coupling 
g^c = VQg°c l or, equivalently, gNNw = 3gNN P - Perhaps a better estimate is obtained 
from the decay width T(uj — > e + e~) = 0.60 ± 0.02 keV, which suggests gNNu> / g NNp = 



[{mJ? p o^ e + e -)/{m p Tu-*e+e-)\ =3.4 ±0.1. 

Keeping only the leading-order contribution proportional to a^ u d u p^, the magnetic cou- 
pling of the vector mesons is defined as f°^/2A4, where the scaling mass Ai, taken to be 



the proton mass, is included to make /° c ct dimensionless. Following Ref. p3| , the SU(6) 



result for a™ can be expressed as a™ = (4M 8 — m l , 8 )/(10M 8 ± 2m v s), where M 8 denotes the 
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average mass in the baryon octet and m„8 the average mass in the vector- meson octet. This 
gives a™ ~ 0.28 for the relativistic SU(6) case, while a™ = | for the static SU(6) case. 

Again applying the VMD hypothesis and assuming that the lowest-mass vector mesons 
(p, d, 0) saturate the nucleon electromagnetic form factors, the magnetic couplings are 
given in terms of the anomalous magnetic moments of the proton and neutron. This gives 
(/ 7 ' 9)nn p = k p — K n — 3.71, and the isoscalar values are expected to be close to (f/g)NNu + 

{f/g)NNcj> ~ Hp + «n = -0.12. 



D. Axial- vector mesons 

We follow Ref. MM in estimating the value for the mixing angle 9 a, defined by 
|/i(1285)) = cos6U K) - sm9 A \a ), 

|/i(1420)> = sm9 A \a$) +cos9 A \a ). (5.15) 

For that purpose we rewrite the /i(1285) and /i(1420) mesons in the nonstrange-strange 
basis, 

| A (1285)) = cos (t) A \A NS ) - sin0 A \A S ), 

|A(1420)> = sin0 A \A NS ) + cos0A \A S ). (5.16) 

The observed decay widths give a mixing angle 4>a ~ 12°. In the singlet-octet basis the 
mixing angle is then given by 9a = 4>a~ arctan \/2 = —42.7°. 

The axial- vector coupling constants are closely related to the pseudovector coupling con- 
stants due to the mixing between the axial-vector and pseudoscalar fields; see the end of 
Sec. |TTT| . The redefinition of the axial-vector fields and the renormalization of the pseu- 
doscalar fields gives the relation between the coupling constants as 

= 9J P- = &9a{0) f 1 -^f] 1 - (5-17) 
The above relation can be recast into the form 



9NNa x — — —fNN-K \ ^ — -, (5.18) 



where the square root equals 1 for Z n = 1/2. The latter choice gives the familiar re- 
sult p3|j35|J45[1 QNNa\ = {nT>ai/ m Tv)fNNTr ~ 8.4. But this assumes m ai = \f2m p , which is not 
supported by experiment. With our previous choice (Z n ,c) = (0.173,-0.131), the axial- 
vector coupling constant comes out to be much larger, w 18. A more moderate value 
is obtained by choosing the alternative solution to Eq. (|3.23| ), (Z^c) = (0.827, 1.25), which 
gives ~ 3.8. An estimate for the a\NN coupling constant from experiment is based 
on the idea of axial-vector meson dominance, which relates it to the axial-vector coupling 
constant of the weak interaction, 5u(0), and the decay constant f ai , as defined by the isovec- 
tor ax-to- vacuum matrix element of the hadronic axial- vector current. With our definition 
( |4.20| ) this gives [|4]] g^f = 4.7 ± 0.6. Hence, we find contradictory results: the larger |c| 



solution gives reasonable agreement with the empirical g^Nai coupling constant, whereas 
the smaller |c| solution gives better phenomenology for a± decay widths and the pion charge 



radius |56]|. This needs to be further explored. 
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VI. DOUBLE-MESON VERTICES 



A. Vector double-meson vertices 



The vector double-meson vertices are obtained from an expansion of Eq. (|4.21b| ). The 
expansion to second order in the meson fields is given by 

$vc = i 7 M (Ap M ) - li \~J 9 g l f f f ^ [(A 77 )' d^XK)] - Y [(Att), (XA,)] + . . . (6.1) 

Hence, it should be obvious that the pair interaction Lagrangians can be obtained by a simple 
replacement of the meson fields in the single-meson interaction Lagrangians. For example, 



the two-pion interaction Lagrangian is obtained from the vector version of Eq. (5A) with 
the replacement 

7 Vm— »7"(txW (6.2) 

which gives 

mlC (n7T) = -g NN ^ n )(NYrN) -(Trxd^) - gsBi^) (E'frE) ■ (tt x <9 m tt) 

-9aa(^) (&fV + S 7 ^A) ■ (tt x d^) + ig^ w) (S x 7 ^S) ■ (tt x d^) . (6.3) 

Here we introduced the square of the charged-pion mass to make the coupling constants 
dimensionless. Substituting the appropriate renormalization factors and \gy = gNNp, the 
coupling constants are given by 

9B'B(nn) = Tf2 ~y ~ > ( 6 ' 4 ) 

4 Ji z tt gNNp 

for B'B = NN, EE, AS, and EE. Note that by choosing |3|j35|j4!| Z w = \, i.e., making 



the assumption that m ai = \[2m p , all the (tttt) pair interactions are absent. This agrees 
with the assumption of vector meson dominance which in its most stringent form states 
that those multiple-meson interactions which can arise through the exchange of one single 
vector meson should not also occur directly P5| . However, experimentally m ai ^ \f2m p , 
and so here the (tttt) pair interactions are still present in the interaction Lagrangian. We 
can also view the (ttit) pair exchange as to proceed via p exchange where the p-meson 
propagator is approximated by a constant (which should be adequate for a heavy meson 
at low energies). Comparing the coupling constant for this and any other pair vertex with 
the effective coupling constant as obtained in such a meson saturation picture then gives an 
indication of how good this picture really is. 

The (ttK) and {rj%K) combinations occur at the same place in the $ vc matrix as the K* 
fields. Therefore, the (itK) and (rj&K) interactions are obtained by substituting, respectively, 

YK; — -iY(Vsd„K - Kdtfa). (6.5) 
This gives the coupling constants 
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— ^ 4V2f? Z K „ g NNp 



m 



2 



2Zk V8 — 1 V~3gB'BK* 



4-\/2/i 9NNp 

for = A7V, SA, EAT, and SE. Here we defined averaged renormalization constants 
Z"^ b = Z a Zb to simplify the expressions. 

In the <3> vc matrix, the two-kaon interactions occur on the diagonal and on the same off- 
diagonal places as the p fields. We can therefore split up the contributions into two parts. 
One part behaves like an isoscalar octet field, and so we have the replacement 

Yuj 8 , — iY{Kd,K - R-^Kc), (6.7) 

and the coupling constants 

ml 2Z K - 1 VSg B Bu;s (a c x 

Qbb(kk) = tth — 7 " • y b - s > 

°Jl gNNp 

The remaining part behaves like an isovector field K 2m which can be written as 

K 2M = (KrdpK) + (K c rd^K c ), (6.9) 

The Lagrangian for this type of two-kaon interaction is obtained by making the substitution 
p^ — > iK 2M , which gives the coupling constants 

ml 2Z K - 1 g B , Bp 

9B'B(KtK) = TTH y ' l b - 1U ) 

°J1 gNNp 

where again B'B = NN, EE, AS, and EE. 

The double-meson interactions consisting of a pseudoscalar and an axial-vector meson 
are analogously obtained by the substitutions 

YPp — ir(KrK lfl + K~ c rK lliC ). (6.11) 

Including averaged renormalization constants, the coupling constants are given by 

m n 

gB>B(ir ai ) — —7=J- gB'Bp, 

m n /y 

gB'B(irKi) — gB' B(Kai) ~ 7=5 , \f ^QB'BK*, 

V ^Kwjl 

m v 

gB'BirisKi) — —gB'B(Kag) ~ 1 — y ^B'BK* ■> 

™* R 

g B' B(K Ki) — „ pj—r v og B >Bujg, 



2VZ^f 
ZVZ^f: 

with the relevant substitutions for the B'B baryon combinations. 



Tfl 

9b'b(KtK\) = 7=7 - gB'Bp, (6-12) 
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B. Axial- vector double- meson vertices 



The expansion of Eq. ( |4.21cj ) to second order in the meson fields gives 

*« = ^7V(A^) + l -Z^^(X.) - ' (1 2 'J / M) 7V [M, (Aft)] + . . . 

(6.13) 

Completely analogous to the previous section, we make the substitutions 

7VA„ — ►tVCtxPaO: 

7 V^ — > -i T V(^ - #u, 8/ J, 
Wo*. — H 5 r(KK; - jQ^J, 

7 5 7^A M — > ^(KtK; + K~ c tK; c ). (6.14) 

The coupling constants are most easily expressed in terms of gmrp and the pseudovector 
coupling constants: 

gB'B(irp) = 2gNNpfB'BTr, 
g B' B(irK*) = —gB'B(Kp) = V2g NNp f B 'BK, 
9B'B(rj a K*) — —9B'B(Klo 8 ) = ^>9NN pf B' BK 
9B 1 B(KK*) = VSgNNpfB'Brjs, 
9B'B(KtK*) — gNN pf B 1 Bit- (6.15) 



C. Possible extensions 



Given the transformation property of the $ sc and $ ax fields, we can of course arbitrarily 
add more interaction Lagrangians of the type ( |4.2Q ) with higher orders in the meson fields. 



In general this will introduce a number of new free parameters. Because there is no guarantee 
that we can actually construct a baryon-baryon potential without any free parameters for the 
coupling constants, it might be very useful to still have some flexibility (i.e., free parameters) 
in the model. It will be convenient, of course, if at the same time we can still satisfy the 
empirical constraints for the single-meson coupling constants as given in Sec. [V|. With this 
purpose in mind, we therefore now discuss some possible extensions in more detail. 

Although most of the scalar meson masses are already close to 1 GeV, the low mass of 
~500 MeV in the two-pole approximation [^9] to the broad £:(760) meson makes the double- 



scalar interactions worthwhile to be investigated. We can extend the scalar interaction 
Lagrangian by adding v / 2(fi , ss/"i7r)( < ^sc) 2 to $ sc of Eq. (|4.21aj) , where the charged-pion mass 



is introduced to make the coupling constant dimensionless. Since the octet and singlet 
parts in the interaction Lagrangian ( |4.20|) have independent coupling constants, we can also 
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define two independent scalar-scalar coupling constants. Hence, writing X = X c s c we have 
the substitution 



S£?*% — ^ (F + X) + H (F 2 + (FX + XF) + X 2 ) | 



(6.16) 



and a similar expression for the singlet part. We still want to satisfy the baryon-mass 
relations of Sec. |V A], which gives the constraints 



flEWi + / 2 ) + **-(2# + £ 



m-n 



(8) 

^°c Ct (/ 2 -/l) + — (/I-/?: 



1153.5 MeV, 
136.5 MeV, 



(6.17) 



where we substituted M = 1153.5 MeV. Proceeding along the lines of the previous section, 
the evaluation of (A c s c ) 2 under the assumption of ideal mixing allows us to read off the 
double-scalar interactions by making the following substitutions 



n oct 
9 SC S 8 



9$ 



2/oao + |(ktk - K C TK C 
(r-a )K + (f - V2e)n 



9^- 



9$ 



4 + ft 



2e 2 - kk 



~ JO 

a o + fo +£ 2 + 2kk 



(6.18) 



where we included the relevant coupling constants. Note that the vacuum expectation matrix 
F does not commute with A c s c , and so also the scalar single-meson coupling constants need 
some modification. In terms of the original coupling constants, we have 



9b'b 



9B'Bh 



9b'b 



*8 



9b'Bs 



9B'Ba 



9B'Bk 



9b'Bs s 



9b'Bs 



(8) 
ss 

oct 



1 + 



m w g 
fi + / 2 9, 



(8) 
.s.s 



rt oct 
i>sc 



(8) 
ss 



1 + 2 fi + 2f 2 g 

3 m n g°™ 
1 + 2 2_h + f 2 g$ 





./s 



2y/2 h-h „ 
9b'Bs s , 



, 2y/2 A - / 2 _ 
H ^ 9b'Bs , 



(6.19) 



where Qb'Bs^ is obtained from the scalar version of Eq. ( |5.3| ) replacing s by sg and gf™ by 
g^l>. Similarly, c/b'Bso is obtained from the scalar version of Eq. (|5.7p replacing s 8 by s 
and g° c ct by gj®. Clearly, this type of extension allows for much more freedom in the scalar 
one-boson exchanges, without the need of abandoning the generation of the proper baryon 
masses. Note, however, that because of the second constraint in Eq. ( |6.17| ), the resulting 
numerical values for gB'Bn, of some importance in YN potentials, remain unaffected by these 
changes. 
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A similar example is the case where we include the matrix for the scalar fields in the 
axial-vector interaction Lagrangian. This gives both scalar-pseudoscalar and scalar-axial- 
vector pair interactions. Since the axial-vector meson masses are well above 1 GeV and the 
scalar meson masses are close to 1 GeV, the scalar-axial-vector exchanges are expected to 
be completely negligible in baryon-baryon potential models, and we will not consider them 
here. The scalar-pseudoscalar interactions are generated by the combination {$ sc ,$ ax }. In 
principle, another possible combination is i[$ sc , $ax], where the i is required by hermiticity. 
However, due to the fact that the vacuum expectation matrix F does not commute with 
$ ax , the commutator generates a complex contribution to the baryon-baryon-kaon coupling 
constants; we therefore drop this combination. As before, by including the anticommutator 
combination we can introduce two coupling constants g^} and g^, for the octet and sin- 
glet part of the interaction Lagrangian, respectively. The form of the various interaction 
Lagrangians and the expressions for the pair coupling constants are now easy to derive. 
Note that the baryon-baryon-pseudovector coupling constants are modified in an analogous 
manner to Eq. Q6.19Q . 

As a final example we discuss the class of pair-meson interactions which, within the 
present formalism, do not have any theoretical constraint on the overall coupling constants. 
The most simple interaction of this type contains ($ ax ) 2 o r ; equivalently, u^u v . This is 
also the most important one since it contains the lightest meson (the pion). Because the 
axial-vector mesons are already rather heavy, in the following we will only consider the 
pseudoscalar-pseudoscalar contributions. We have the possibility for two types of field com- 
binations, one symmetric and one antisymmetric in the fields: 

<p s ~ -§<T [0 M (A7r)0„(A7r) + ^(A^Atf)] , 

0a ~ +f ^ [^(AtO^Att) - d„(\n)d^\n)] . (6.20) 

It will be convenient to identify the two-pseudoscalar contributions with a matrix of scalar 
fields as in Eq. (|4.14|) . For the symmetric combination this implies the following substitutions 



a - 




K - 




S8 ~ 




so - 





\t- [d^dyK + dvirdpK) - J\{d^d v K + d u f] 8 d^K) + Jj(d fl r} d v K + d^d^K) 



2 

\d^Tz-d u Tz - J\d^r] 8 d v r]s 



dfx'K • d u 7r + df.rjsd^s + d^d^o + (d^Kd u K + d u Kd^K) 
while for the antisymmetric combinations we have 

id^n x d u ir - \(d,JCTd v K - dyKrd^K) 



(6.21) 



a 

K 
S8 



\r-{d^d v K - d u ird^K) 
^{d^K-dJKd^K) 



(6.22) 



These matrices can be substituted in the interaction Lagrangian ( |4.20| ), which contains the 
free parameters gj$ m , g$ m , and a sym for the symmetric case, and g£j ym and a asym for the 
antisymmetric case. Including the renormalization factors it is now straightforward to find 
the pair coupling constants for each of the two-pseudoscalar contributions, expressed in 
terms of these free parameters. 
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VII. APPLICATION TO NN 



As a first application of the chiral-symmetry constraints given in this paper, we like 
to investigate whether with the values for the coupling constants as given in the previous 
sections it is indeed possible to construct a baryon-baryon potential model which gives 
a satisfactory description of the baryon-baryon scattering data. Of course, the imposed 
constraints need not all be exactly true. For example, the vector- dominance assumption 
that k p = 3.71 and + — —0.12 is only true if these mesons fully saturate the nucleon 
electromagnetic form factors. The presence of heavier vector-meson nonets likely changes 
these relations P5 . Also, the SU(3) relations in Sec. [V] need not be true in an exact sense. 
This is already clear from the fact that we have to introduce symmetry-violating terms 
to generate the empirical meson masses. The existence of a scalar meson nonet and its 
quark content is still under debate, and so the assumption of an SU(3) symmetry for the 
scalar mesons might even be incorrect. On the other hand, relaxing too many constraints 
introduces too many free parameters - something we would like to avoid. Therefore, at this 
stage we choose to impose all the constraints and only show that the resulting NN potential 
model then already gives a very reasonable description of the scattering data. 

The experience with NN potential models that have appeared in the literature sug- 
gests that a fully constrained potential model of the one-boson-exchange type is unlikely 
to succeed [|]. On the other hand, we have already demonstrated pT|JT^] that by includ- 
ing two- meson-exchange contributions a major improvement in the description of the NN 
scattering data can be obtained. In order to arrive at a model which at least gives a reason- 
able description of the scattering data, we found that we had to include the double-scalar 
and double-pseudoscalar extensions as outlined in Sec. |V1 C[ The single-meson coupling 
constants satisfy the empirical constraints as discussed in Sec. [V| and the pair-meson cou- 
pling constants satisfy the relations as given in Sec. |VT|. The one-boson-exchange part of 
the potential is standard but includes the diffractive contribution |[47|| , while the two-meson 
part can be found in, or easily derived from, Refs. |Iiyi2"fl . The potential is regularized 
with exponential form factors, one for each type of meson (scalar, pseudoscalar, or vector). 
The single-meson coupling constants and the exponential cutoffs for each type are given in 
Table |. The pair-meson coupling constants are given in Table ||. Note that we only include 
meson pairs with a total mass below ~1 GeV. However, since the ^-exchange contributions 
did not significantly improve the fit, we decided not to include them at this stage. 

The 12 free parameters of the model (A 5 , A P , Ay, #a 2 , 9p, gj$, g&\ gi%, a sym, 
fidsj/m.5 an d OLasym) were determined in a fit to the Nijmegen representation jPJ of the \ 2 
hypersurface of the NN scattering data below T lab = 350 MeV, updated with the inclusion 
of new data which have been published since then. The effective diffractive mass, mp = 310 
MeV, was fixed at the (rounded-off) value as used in the old Nijm78 potential ^7|. The 
resulting xV^data for each of the ten energy bins is shown in Table |T|, in comparison with 
the (updated) Nijmegen partial-wave analysis. 

The xV^data = 1-75 for the 0-350 MeV energy interval actually compares very favorably 
to other potential models that have appeared in the literature 0. As a matter of fact, it 
should be realized that in this model all coupling constants satisfy constraints as imposed 
by chiral symmetry, or empirical constraints as discussed in Sec. |V|. This in contrast to 
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any other model that has appeared in the literature. The model even gives a much better 
description of the data below 300 MeV (xV^data = 1-36), whereas it rapidly worsens at 
higher energies. This sudden rise is probably due to the nonadiabatic expansion in the two- 
meson contributions [fTT| , [T2|1 which, strictly speaking, is only valid below the pion production 
threshold (Tj a b ~ 280 MeV). The nonadiabatic expansion is an artifact of us working in 
coordinate space, and the sudden rise in x 2 a t higher energies will be further investigated 
when we have developed a momentum-space version where we can retain the full energy 
dependence in the propagators^. At this stage we prefer to work in coordinate space for 
several reasons. First, the (already rather time consuming) fit in coordinate space is much 
faster than a fit in momentum space. Second, here we only wanted to investigate whether 
it is indeed possible to construct a potential model which incorporates all constraints and 
still gives a satisfactory description of the data. Modifications such as keeping the full 
energy dependence in the propagators rather than making a nonadiabatic expansion can 
be investigated at a later stage. Third, ultimately it is our goal to produce a high-quality 
potential model which is exactly equivalent in both coordinate space and momentum space. 
This allows equivalent applications in both spaces, but requires certain approximations to 
be made in order to arrive at analytical expressions in coordinate space. 

We should point out that, although the number of free parameters is not too large, the 
fully constrained fit is far from trivial. Due to the constraints, the relation between the 
change in a parameter and the corresponding change in the phase shifts is highly nonlinear. 
It turned out to be impossible to fit all the free parameters at the same time, and so we 
had to do numerous fit cycles where we fit only an (arbitrary) subset of the parameters. 
This makes it very complicated and time consuming (but not impossible) to arrive at a 
satisfactory fit, and we cannot guarantee that the fit presented here is the most optimal 
one. However, the present result already clearly illustrates our main objective, viz. to try 
to construct a potential model which gives a reasonable description of the scattering data, 
and at the same time incorporates a number of empirical and chiral-symmetry constraints. 
Improvements by investigating different parameter sets, adding more parameters, or relaxing 
some of the constraints will be left for the future. 

VIII. SUMMARY 

We have constructed a chiral-invariant Lagrangian for the meson-baryon sector, where 
the mesons consist of the nonets of scalar, pseudoscalar, vector, and axial-vector mesons, 
and the baryons are the members of the baryon octet. Although we mainly focussed on 
the meson-baryon interaction Lagrangian, we briefly indicated what the Lagrangian in the 
meson sector looks like, and how it can be made to reproduce the correct phenomenology 
by allowing for some small contributions which violate the chiral symmetry. 



Alternatively, we can decide that the model should only be used up to T^b ~ 300 MeV, say. 
Without pursuing this option any further, a quick and not too thorough refit already shows that 
a x 2 /Vdata = 1-3 seems easily feasible. 
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In the meson-baryon sector, the chiral symmetry imposes constraints on the coupling 
constants of the meson-pair vertices. This allows us to express all meson-pair coupling 
constants in terms of the single-meson vertex coupling constants. These single-meson vertex 
coupling constants, in principle, can all be fixed by using experimental data such as meson 
masses, baryon masses, meson decay parameters, and meson mixing angles. 

As a first step in the application of the meson-baryon Lagrangian, we demonstrate that 
it is possible to construct an AW potential which has all the theoretical and empirical 
constraints as discussed in Sees. [V| and [Vl|. The model gives a very satisfactory description 
of the AA scattering data. It even has a slightly better quality than other (often more 
phenomenological) potential models that have appeared in the literature, where this is the 
first model where such strict constraints on the coupling constants have successfully been 
imposed. We should mention that this result could only be achieved after the inclusion of 
the two-meson contributions. Improvements and an extension to the YN sector are presently 
under investigation. The success of the present model gives us hope that it might indeed 
be possible to ultimately arrive at a potential model which not only provides a high-quality 
description of the scattering data, but which is also consistent with the symmetries of QCD. 
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APPENDIX: 



We follow Ref. J|9| and show that it is possible to find a nonlinear transformation that 
transforms away the octet of original pseudoscalar fields, while leaving the scalar vacuum 
expectation matrix F invariant. We write the scalar nonet as X c a c = F + E, where the 
diagonal matrix F = (fi, f'2, ^3) is the nonvanishing vacuum expectation matrix of the scalar 
fields. The parameters /, can be anything, but in the following we assume that fi + f) 7^ 
for all {i, j}. The pseudoscalar octet is written as II = A a 7r a . We want to find a traceless 
matrix P (to be identified as an octet of new pseudoscalar fields), such that 

F + £ + iU = e ip (F + X + iY)e iP } (Al) 

where X and Y are sets of different scalar and pseudoscalar fields, respectively. For that 



purpose, we formally write [[49 
P 

and try to solve for 



00 

a a k P k 

k=0 



X 



k=0 



Y 



k=0 



(A2) 



F/a + E + iU = e iP (F/a + X + iY)e iP , (A3) 
in each order of a. To order a -1 we have an identity. To order a we can write 
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v _ v v _ nil/2/3 + II22/1/3 + ^33flf2 fl 
X0 " S ' / 2 /3 + /l/3 + /l/2 

UU ; (/< • /,) ' II Vo, r (A4) 

where we chose Yq such that Tr(P) = 0. This is not the only possible choice for Yq, but it 
will turn out to be a very convenient one. (Note that in Ref. |4£| a different choice is made, 
but there Pq is not traceless, contrary to what is claimed.) This procedure can be extended 
to all orders in a, and we always find that 

Y n + {F,P n } = g n , (A5) 

with g n a nonlinear function of Pf., X^, and Y^ with k = 0, . . . , n — 1. Hence, in analogy to 
the n = case, we can always define a traceless matrix P n and at the same time keep Y n 
proportional to H 3 . Substituting u = exp(zP), the definition (|4.8Q gives 



X+ = F + X, X -=^, (A6) 

where X represents the new scalar octet and where Y (because it is proportional to H 3 ) can 
be identified with a new isosinglet pseudoscalar field, not present before. Clearly, different 
choices for Y n , not proportional to I3, will introduce unwanted contributions to the isovector 
(71*) and/or isoscalar (77s) pseudoscalar fields. 
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TABLES 



TABLE I. Single-meson NN coupling constants and exponential cutoff masses for the con- 
strained NN model. The single-meson coupling constants are divided by ^/4tt. 





g-scalar 


/-pseudosc. 




g- vector 


/-vector 


g-diffr. 




£ /o 


tt rj rj' 


P 


LO <j) 


p UJ <p 


A 2 P 


2.599 


3.273 -0.978 


0.270 0.117 0.075 


0.711 


2.408 -0.183 


2.634 0.297 0.045 


2.497 2.126 


A 5 : 


= 590.0 MeV 


A P = 870.0 MeV 


Av 


= 845.0 MeV 







TABLE II. Pair-meson coupling constants for the constrained NN model. The coupling con- 
stants are divided by An. For the double-pseudoscalar irir and KK contributions we distinguish 
interactions symmetric and antisymmetric in the fields. 













symmetric 


antisymmetric 


7T7T 


KK KtK irp 


EE 


177] 


ITTj' 


7T7T KK KtK 


7T7T KK KtK 



-0.030 -0.028 -0.009 0.384 2.463 0.007 0.004 0.045 0.007 0.005 -0.050 0.082 0.033 



TABLE III. x 2 an d X 2 P er datum (x 2 , d p ) a ^ ^ ne 10 energy bins for the updated partial-wave 
analysis (PWA) and the constrained NN potential. A^ata lists the number of data within each 
energy bin. The bottom line gives the results for the total 0-350 MeV interval. 



PWA potential 



Bin(MeV) 


^data 


x 2 


Xp.d.p. 


x 2 


Xp.d.p. 


0.0-0.5 


145 


144.45 


0.996 


155.73 


1.07 


0.5-2 


68 


43.08 


0.633 


56.48 


0.83 


2-8 


110 


105.11 


0.956 


157.41 


1.43 


8-17 


294 


275.24 


0.936 


335.07 


1.14 


17-35 


359 


294.58 


0.821 


387.39 


1.08 


35-75 


585 


565.30 


0.966 


1045.22 


1.79 


75-125 


399 


409.62 


1.027 


443.18 


1.11 


125-183 


760 


823.82 


1.084 


1249.44 


1.64 


183-290 


1046 


1022.72 


0.978 


1306.90 


1.25 


290-350 


992 


993.49 


1.002 


3197.26 


3.22 


0-350 


4758 


4677.42 


0.983 


8334.09 


1.75 
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